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ABSTRACT 

IRAS 19475+3119 is a quadrupolar pre-planetary nebula (PPN), with two 
bipolar lobes, one in the east-west (E-W) direction and one in the southeast- 
northwest (SE-NW) direction. We have observed it in CO J=2-l with the Sub- 
millimeter Array at ~ 1" resolution. The E-W bipolar lobe is known to trace 
a bipolar outflow and it is detected at high velocity. The SE-NW bipolar lobe 
appears at low velocity, and could trace a bipolar outflow moving in the plane 
of the sky. Two compact clumps are seen at low velocity around the common 
waist of the two bipolar lobes, spatially coincident with the two emission peaks 
in the NIR, tracing dense envelope material. They are found to trace the two 
limb-brightened edges of a slowly expanding torus-like circumstellar envelope 
produced in the late AGB phase. This torus-like envelope originally could be 
either a torus or a spherical shell, and it appears as it is now because of the two 
pairs of cavities along the two bipolar lobes. Thus, the envelope appears to be 
disrupted by the two bipolar outflows in the PPN phase. 

Subject headings: circumstellar matter - planetary nebulae: general - stars: AGB 
and post-AGB - stars: individual (IRAS 19475+3119): mass loss -stars:winds, 
outflows 



1. Introduction 



Pre-planetary nebulae (PPNe) are objects in a transient phase between the asymptotic 
giant branch (AGB) phase and the planetary nebula (PN) phase in the stellar evolution of 
a sun-like star. The PPN phase is short with a lifetime of up to thousands of years; yet it 
closes the gap between the two morphologically very different phases - the AGB phase and 
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the PN phase. The former exhibits roughly spherical AGB circumstellar envelopes (CSEs) 
expanding radially at 10-20 km s _1 , comparable to the escape velocity of the AGB stars. Yet 
the latter exhibits diverse morphologies, with a significant fraction having h ighly collimated 
fast-moving bipol ar or multipolar (including quadrupolar) outflow lobes (IBalick fc Frank 
20021 ; ISahail I2004J ) . In the past decades, high spatial resolution imagings of PPNe e xhibit 
an even richer array of shapes, including multiaxial symmetries ( IBalick &: Frankj|2002l ). sug- 
gesting that the aspherical symmetry in the PN phase must have started in the PPN phase 
( ISahail 120041 ) . It seems that collimated fast winds or jets are required to shape the PPNe 



and young PNe (Sahai fc Traugerlll998l ; ISoker fc Rappaportll2000l ; iLee fc Sahaill2003t ISoker 
2004 ; ISahail 120041 ) . In particular, the departure from the spherical symmetry could result 
from an interaction of the jets with the preexisting AGB CSEs. Also, multipolar and/or 
point-symmetric morphologies could result from multiple ejecti ons and /o r change in ejection 
directions of the jets, respectively (Sahai et al. 2007, hereafter ISSMC07I ) . 



PPNe and PNe are often seen with a dense equatorial waist structure that is ofte n 
called a torus (e.g., 



Kwok et al. 


1998; 


Su et al. 


1998; 


Volk et al. 


2007; 


Sahai et al. 


2008) 



The origin of the torus is unknown, and it may be due to either an equatorially enhanced 
mass loss in the late AGB phase or an interaction between an underlying jet and the spherical 
AGB wind ( ISoker fc Rappaportll2000l ). Tori are common in bipolar PPNe, but less common 
in multipolar PPNe. It is thus important to study how a torus can be formed in multipolar 
PPNe and its relation with the jet. 

The PPN IRAS 19475+3119 is a quadrupolar nebula with two collimated bipolar lobes, 
one i n the east -west (E-W) direction and one in the southeast-northwest (SE-NW) direc- 
tion (ISSMC071 ). In the optical image of the Hubble Space Telescope (HST), the bipolar 
lobes appear limb-br ightened, suggesting that they are dense- walled structures with tenuous 
interiors or cavities ( 1SSMC07I ). and that they trace two bipolar outflows produce d by two 
bipolar post-AGB winds (Sanchez-Contreras et al. 2006, hereafter ISC06I ; ISSMC07I ). Its cen- 
tral star, HD 331319, is an F3 lb (T e ff ~ 7200K) s upergiant star classified as a post-AGB 
object based on the elemental abundance analysis (jKlochkova et all 120021 ). The distance 
is ~ 4.9 kpc by assuming a total luminosity of 830 L , appropriate for a post-AGB star 
with a mass of 0.63 M ( IHrivnak fc Biegingl 120051 ). Using this luminosity and the effec- 
tive temperature, the central star is estimated to have a radius of 58 R . This PPN has 
a detached circumst ellar envelope based on its double-peaked spectral energy distribution 
( IHrivnak et al.l 119991 ) . Near-infrared (NIR) images of this PPN show a nebula extending i n 
the SE-NW direction and a puzzling two-armed spiral-like structure (IGledhill et al.l 120011 ). 
The NIR nebula is the counterpart of the SE-NW bipolar lobe in the optical. The spiral-like 
structure is point-symmet ric and could result from an interaction of the mass- losing star 
with a binary companion (IGledhill et al.ll200ll ). 



-3 - 



Single-dish CO spectra toward this PPN show two components, a strong line core 



from the circumstellar e nvelope and weak wings from fast moving gas (ILikkel et al.l 11991 
Hrivnak fc Biegingl 120051 ) . At ~ 2" resolution, the two components are seen in CO J=2-l 



as a slowly expandi ng she ll and a fast bipolar outflow that is aligned with the E-W bipolar 
lobe in the optical (1SC06I ). At ~ 0"6 resolution, a ring-like en velope is seen mainly around 
the E-W bipolar lobe ( Castro- Carrizo et al. 2010, hereafter ICClOl ). Here we present an 
observation at ~ 1" resolution in CO J=2-l for this PPN. Our observation shows an ex- 
panding torus-like envelope around the common waist of the two bipolar lobes and that the 
structure and kinematics of the envelope are consistent with the envelope being disrupted 
by the quadrupolar nebula, or the two bipolar outflows. 



2. Observations 

Interferometric observati on of IRAS 194 75+3119 in CO J=2-l was carried out using the 



Submillimeter Array (SMA) ( )Ho et al.l 120041 ) in the extended configuration on 2009 August 
29. Seven antennas were used with baselines ranging from 44 to 226 m. The SMA correlators 
were set up to have a channel spacing of 0.2 MHz (~ 0.26 km s _1 ) for CO J=2-l. The total 
duration of the observation (including the integration time on the source and calibrators) 
was ~ 11.5 hr with a useful duration of ~ 5.5 hr on the source. 

The calibration of the data was performed using the MIR software package. Passband 
and flux calibrators were 3C84 and Uranus, respectively. The gain calibration was done with 
the quasar 2015+371, whose flux was estimated to be 3.05 Jy at 230 GHz, about 20% higher 
than that reported from the light curve measurements of the SMA around the same period of 
time. Each integration cycle includes 4 minutes on 2015+371 and 20 minutes on the source. 
The calibrated visibility data were imaged with the MIRIAD package. The dirty maps that 
were produced from the calibrated visibility data were CLEANed, producing the CLEAN 
component maps. The final maps were obtained by restoring the CLEAN component maps 
with a synthesized (Gaussian) beam fitted to the main lobe of the dirty beam. With natural 
weighting, the synthesized beam has a size of l"16xl // 01 at a position angle (P. A.) of —63°. 
The rms noise levels are ~ 0.1 Jy beam" 1 in the CO channel maps with a channel width of 
1 km s" 1 . The velocities of the channel maps are LSR. 
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Results 



In the following, o ur results are p r esente d in comparison to the HST image of ISSMC07 



and the NIR image of iGledhill et all ( 1200 If ). Note that the HST image adopted in ISC06 
was incorrectly rotated by ~ 10° clockwise. The systemic velo city in this PPN is found to 



be 18.25±1 km s" 1 from our modeling in Section H] (see also ISC06I ). The integrated line 
profile of CO shows that most emission comes from within ±15 km s _1 from the systemic 
velocity. Com pared with previous single-dish and interferometric observations in the same 
emission line ( 1SC06I ). our interferometric observation at high resolution filters out most of 
the extended emission at low velocity, allowing us to study compact structures, such as a 
torus and collimated outflows, around the source. 



3.1. A torus-like circumstellar envelope 

At low velocity (within ± 10 km s _1 from the systemic velocity), the CO map (Fig. [I]) 
shows two compact clumps on the opposite sides of the source, one to the northeast (NE) 
at PA. ~ 15° and one to the southwest (SW) at PA. ~ 195°, at a distance of ~ 1" (4900 
AU) from the source. These two clumps can be seen in the 7 central ve locity channels from 



8 to 29 km s _1 (Fig. [2]) and have also been seen at higher resolution in lCClOl as well. They 
are seen around the common waist of the two bipolar lobes, spatially coincident with the 
two emission peaks in the NIR, tracing the dense envelope around the source. The center 
of the two clumps, however, is shifted by Of' 15 to the east of the repo rted source position, 



which is a (20 oo) = 19 h 49 m 29 s 56, 5 (2 ooo) = +31°27'16"22 f lHog et al.lll998f ). This position shift 
is only ~ 15% of the synthesized beam and thus can be due to the position uncertainty in 
our map. Also, there could be uncertainty in the reported source position. Therefore, in this 
paper, the center of the two clumps is considered as the source position. Faint protrusions 
are also seen along the two bipolar lobes at low velocity, extending to the east and west, as 
well as to the southeast and northwest. These faint protrusions are most likely to be from 
the quadrupolar lobes, and thus unlikely to be from the dense envelope itself. 

The structure and kinematics of the two clumps can be studied with position-velocity 
(PV) diagrams cut along four different PA., from the one aligned with their peaks (cut A 
with PA. =15°) to the one perpendicular (cut D with PA. = —75° ), through the two in 
between (cut B with PA. = — 10° and cut C with PA. =—30° ). The PV diagrams along these 
cuts shows a similar ring-like PV structure with less and less emission around the systemic 
velocity as we go from cut A to cut D (Figures Ek-d), indicating that the two clumps trace 
the two limb-brightened edges of a torus-like envelope at a small inclination angle (i.e, close 
to edge-on), either radially expanding or collapsing, with the equator aligned with cut A. In 
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addition, the torus-like envelope has a thickness with a half opening angle (i.e., subtended 
angle or see next section for the definition) of 20°-45°, measured from its equator. The 
two high-velocity ends of the ring-like PV structure are from the farside and nearside of the 
torus-like envelope projected to the source position. The emission there is bright, indicating 
that the material in the farside and nearside must be projected to the source position, and 
this requires the torus-like envelope to have an inclination angle smaller than its half opening 
angle. Figure H] shows the spectrum integrated over a region with a diameter of ~ 2"5 around 
the source. It includes most of the torus-like envelope but excludes the protrusions along 
the two bipolar lobes. It shows a double-peaked line profile with the redshifted peak at ~ 28 
km s _1 slightly brighter than the blueshifted peak at ~ 10 km s _1 . Since t he temperature 
of th e torus-like envelope is expected to drop away from the central star ( IK wan fc Linke 



19821 ). the emission from the farside is expected to be brighter than that from the nearside 



because of the absorption in the nearside. Therefore, the redshifted peak in the line profile is 
from the farside and the blueshifted peak is from the nearside, implying that the torus-like 
envelope is expanding at ~ 10 km s" 1 , in agreement with it being the circumstellar envelope 
produced in the AGB phase. The PV diagram cut perpendicular to the equator, i.e., cut 
along P.A.=-75°, shows that the redshifted emission shifts to the west and the blueshifted 
emission shifts to the east at low velocity (Fig. [3U). For an expanding torus, this indicates 
that the torus is inclined, with the nearside tilted to the east. 



3.2. Outflows 



At high velocity (beyond ±13 km s~ l from the systemic velocity), CO emission is seen 
around the E-W bipolar lobe seen in the optical, emerging from the two opposite poles (or 
holes) of the torus-like envelop e (Fig. [T]). This E-W bipolar lobe is already known to trace 
a fast-moving bipolar outflow (jSC06l ). The PV diagram (Fig. |3fe) along the outflow axis 
(P.A.=88°; cut E in Fig. [1^) shows th at the outflow velocity increases from the base to the 
tip at ~ 4" (~ 20000 AU), as found in lSC06l . This is why the CO emission at low- velocity is 
seen near the source forming the faint protrusions along the E-W bipolar lobe, as mentioned 
in the previous section. The radial velocities are ~ 30 km s -1 at the tip of the red-shifted 
lobe and ~ 24 km s _1 at the tip of the blue-shifted lobe, with a mean velocity of ~ 27 km 



The PV diagram also shows a hint of bifurcation in velocity along the outflow axis, 
especi ally fo r the red-shifted lobe in the east [see the dark lines in Figure |3k and also Figure 



20 in |CC10I |. This su ggests that the outflow lobe is a hollow, shell-like cavity wall, as 



suggested in ISSMC07I . For example, in the case of the redshifted lobe, the emission with 
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less red-shifted velocity is from the front wall, while the emission with the more red-shifted 
velocity is from the back wall. The outflow has the shell-like structure, likely resulting from 
an intera ctio n of a pos t-AGB wind with the pre-existing circumstellar medium, as proposed 



in 



SC06 and SSMC07. 



On the other hand, no CO emission is seen here at high velocity around the SE-NW 
optical bipolar lobe. It has been argued that this bipolar lobe, lik e the E-W bipolar lobe, is 
also a dense-walled structure with tenuous interiors (or cav ities) ( 1SSMC071 ) and could also 
be produced by a bipolar post-AGB wind ( ]SC06l ; ISSMC071 ). Thus, this bipolar lobe could 
trace a bipolar outflow moving along the plane of the sky, with the CO emission seen only 
at low velocity (within ±10 km s _1 from the systemic as seen in the channel maps) forming 
the faint protrusions along the SE-NW bipolar lobe, as mentioned in the previous section. 
The low inclination is also consistent with the fact that the SE and NW components of the 
SE-NW bipolar lobe have almost the same brightness in the optical. 



4. Modeling the torus-like envelope 

As discussed, the two CO clumps seen at low velocity arise from a torus-like envelope 
at a small inclination angle, with the nearside tilted to the east. In this section, we derive 
the physical properties, including the kinematics, structure, density and temperature distri- 
butions, of this torus-like envelope through modeling the two CO clumps with a radiative 
transfer code. The envelope appears torus-like around the common waist of the two bipolar 
lobes, and it could arise from a torus with or without the cavities (or holes) cleared by the 
bipolar lobes, or a spherical shell with the cavities. For a torus and a spherical shell, the 
number density of the molecular hydrogen in spherical coordinates (r, 9, 0) can be assumed 
to be _ 2 

n = n J— ) cos p 6, (1) 
\r J 

where 9 is measured from the equatorial plane of the envelope. Here p = for a spherical 
shell and p > for a torus with a half opening angle 9q defined as 

cos p # = 0.5 (2) 

Also r =1", which is the representative radius of the torus-like envelope. The envelope has 
an inner radius R\ n and an outer radius -R ou t- It is expanding radially at a velocity of v e , and 
thus the mass-loss rate (including helium) can be given by 

M = 1.4m H2 J nv e r 2 d^l = 2.87rm} l2 n v e rl J cos p 9 d9 (3) 
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The cavities, if needed, are assumed to have a half-ellipsoidal (like paraboloidal) opening for 
simplicity, with the semi-major axis and semi-minor axis determined from the bipolar lobes 
in the optical. Inside the cavities, the number density of the envelope is set to zero. 



In the envelope, the temperature profile is assumed to be 

t = tJ- 

Vr 



similar to that derived by iKwan fc Linkd (119821 ) for AGB envelopes, including 
collisional heating, adiabatic cooling, and molecular cooling. 



(4) 

gas- dust 



In our models, radiative transfer is used to calculate the CO emission, with an assump- 
tion of local t herma l equilibrium. The relative abundance of CO to H2 is assumed to be 
4xl0~ 4 , as in |SC06| . For simplicity, the line width is assumed to be given by the thermal 
line width only. Also, the systemic velocity v sys is assumed to be a free parameter. The 
channel maps of the emission derived from the models are used to obtain the integrated 
intensity map, spectrum, and PV diagrams to be compared with the observation. Note that 
the observed emission in the east, west, southeast, and northwest protrusions at low velocity 
(Fig. [T]) is not from the envelope and thus will not be modeled here. In our models, the free 
parameters are v sys , v e , R im R ut, P (torus or spherical shell), cavities (holes), inclination 
angle i, equator PA, n , and T . Table [1] shows the best-fit parameters for our models with 
different structures. Our models all require v sya ~ 18.25±1 km s _1 and v e ~ 12.5±1.5 km 
s _1 to fit the two velocity ends in the PV diagrams and the spectrum. They require Rj n ~ 
0"7 (3430 AU) and R out ~ lf'6 (7840 AU) to match the emission peak positions. They also 
require T to be ~ 28 K to match the emission intensity, resulting in an overall character- 
istic (or mean) temperature of ~ 23 K for the envelope. This character istic t emp erature is 
slightly higher than the values found in the more extended envelope by ISC06I and ISSMC07 
at lower resolutions. This is reasonable, because the envelope temperature is expected to be 
higher near the source. In the following, we discuss our different models in detail. 



4.1. A Torus without cavities 

First we check if a simple torus model can reproduce the two clumps (Figure [5]). In 
this model, the equator of the torus is set to be at P.A.=15°, to be aligned with the two 
clumps. As expected, the torus is required to have a small inclination angle of ~ 15° and a 
small half opening angle of 6$ ~ 24° with p = 8. This model can reproduce reasonably well 
the compactness of the two clumps. It can also reproduce the ring-like PV structures, and 
the less and less emission around the systemic velocity in the PV diagrams from cut A to 
cut D. Note that the PV diagram cut along the equatorial plane shows a ring-like structure 
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with two C-shaped peaks at the two high velocity ends facing to each other, because the 
emission is symmetric about the equatorial plane. For the PV diagrams with a cut away 
from the equatorial plane, the blueshifted emission shifts to the southeast and the redshifted 
emission shifts to the northwest, because the torus is inclined with the nearside tilted to the 
southeast. For the spectrum, the model can produce the double-peaked line profile, however, 
it produces higher intensity than the observed around the systemic velocity (Fig. |5b). 



4.2. A Torus with a E-W pair of cavities 

Here we add a E-W pair of cavities to the above simple torus model and see if we 
can reproduce the spectrum better by reducing the intensity around the systemic velocity 
(Figure [6]). The E-W pair of cavities are assumed to have a half-ellipsoidal opening with a 
semi-major axis of 2" and a semi-minor axis of 1", as measured from the HST image (Fig. [1]). 
The E-W pair of cavities have a P.A.=88° and can have an incl ination angle from 10° to 30°, 



with the western cavity titled toward us, as suggested in ISC06I . These moderate inclinations 
are consistent with the optical image that shows the eastern lobe is slightly fainter than 
the western lobe, and they all give similar results in our model. The model results here are 
similar to those of the simple torus model above, with some slight differences. This is because 
the E-W pair of cavities are almost aligned with the poles of the torus, and removing the 
tenuous material in the poles does not change much the model results. The cavities remove 
slightly the blueshifted emission in the west and redshifted emission in the east, mainly at 
low velocity, as seen in the spectrum. In addition, unlike that in the torus model, the PV 
diagram cut along the equatorial plane shows a ring-like structure with two sickle-shaped 
peaks, instead of two C-shaped peaks, at the two high velocity ends facing to each other 
(comparing Figure \&p and Figure Et). This is because that, along that cut direction, the 
E-W pair of cavities remove the redshifted material in the north, and blueshifted emission 
in the south. 



4.3. A Torus with 2 pairs of cavities 

Here we add one more pair of cavities, the SE-NW pair, into our model (Figure [7]). 
Again, this pair of cavities are assumed to have a half-ellipsoidal opening but with a semi- 
major axis of 2" and a semi-minor axis of Of'8, as measured from the HST image (Fig. [I]). This 
pair of cavities have a P.A.=— 40° and are assumed to have an inclination of 0° (i.e., in the 
plane of the sky). Adding this pair of cavities will rotate the two clumps counterclockwise 
about the source. Therefore, the torus in this model is required to have its equator at a 
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smaller P.A. of 10° in order for the two clumps to appear at P.A.=15°. As a result, the 
underlying torus could be initially more perpendicular to the E-W bipolar lobe. Also, a 
higher density is required to reproduce the same amount of flux in the clumps (see Tabled]). 

By removing the material along the SE-NW bipolar lobe, this model can now reproduce 
the required dip around the systemic velocity in the double-peaked spectrum. This is because 
that the SE-NW pair of cavities, being in the plane of the sky, remove the low-velocity 
emission preferentially. This model can also reproduce the ring-like PV structures, and 
the less and less emission around the systemic velocity in the PV diagrams from cut A 
to cut D. Note that, in the PV diagram cut along the equatorial plane, the model shows 
two sickle-shaped peaks, which could be somewhat different from the observation in detail. 
This detailed difference could be due to localized excitation effect near the base of the E-W 
cavities and should not affect our main conclusions on the envelope properties. 

4.4. A Spherical shell with 2 pairs of cavities 

A spherical shell with only the E-W pair of cavities can be ruled out because it would 
produce two clumps that are exactly perpendicular to the E-W bipolar lobe, inconsistent 
with the observation. Therefore, here we check if a spherical shell with the same 2 pairs of 
cavities as above can produce the two clumps in the observation. We find that this model is 
quite similar to the torus model with the 2 pairs of cavities, except that the two clumps are 
slightly rotated and the emission is more extended perpendicular to the torus-like structure, 
in between the cavities (comparing Figure Eh to Figure [T^l). This slight rotation is acceptable 
because our model is simple and the two clumps are not resolved in our observation. On the 
other hand, further observation that can separate the envelope from the two bipolar lobes 
is needed to check if the emission extending perpendicular to the torus-like structure in this 
model is inconsistent with the observation. 



4.5. Model summary 

Our model results show that the two CO clumps could arise from either a torus or a 
spherical shell, with two pairs of cavities, one along the E-W bipolar lobe and one along 
the SE-NW bipolar lobe. A pure torus model and a torus with only a E-W pair of cavities 
both produce more than observed emission around the systemic velocity in the spectrum. In 
addition, the SE-NW pair of cavities, being in the plane of the sky, are essential to remove 
the low-velocity emission in our models. All the models are optically thin. With the mean 
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envelope temperature of 23 K and the model brightness temperature of ~ 4 K, the mean 
optical depth is ~ 0.2. Therefore, in the model spectra, the redshifted peak is only slightly 
brighter than the blueshifted peak, as observed. The mass-loss rate of the AGB wind derived 
from our two best models is ~ 2 — 3 x 10~ 5 M yr" 1 . In the shell model, the mass-lo ss rate 
is about three times smaller than that derived from a similar shell model by ISC06I . which 
was found to be ~ 10 -4 M yr" 1 . This is probably because of the following reasons: (1) 
our inner radius is smaller than theirs, and the emission is more efficient in the inner region 
where the density and temperature are both higher; (2) a large amount of missing flux in 
our observation; and (3) their envelope emission seen at low resolution could be significantly 
contaminated by the outflows. Therefore, the mass loss rate derived here could be a lower 
limit. Alternatively, the mass-loss rate might have decreased with time. 



Discussion 



5.1. The expanding torus-like envelope 



In our observation, two clumps are seen in CO at the common waist of the two optical 
bipolar lobes, spatially coincident with the two peaks in the NIR polarization image (Figure 
[J), tracin g the dense envelope material. These two clumps are a lso seen at higher resolution 
by ICClOl . but not spatially resolved at lower resolution in |SC06| . From our models, they are 
found to trace the two limb-brightened edges of an expanding torus-like envelope around 
the source, with an inclination angle of ~ 15°. This torus-like envelope corres ponds to the 
ring-like envelope found at higher resolution at a similar inclination angle in ICClOl . At a 
lower resolution of ~ 2", tw o different clumps were seen in CO but aligned with the SE-NW 
optical bipolar lobe ( 1SC06I ) as well as the SE-NW NIR nebula. Those clumps are also seen 
here but appear as the two faint protrusions along the SE-NW bipolar lobe (Fig{Th). As 
mentioned before, the two optical bipolar lobes, both E-W and SE-NW, h ave been found to 
be dense- walled structures with tenuous interiors (i.e., cavities) (ISSMC071 ) and these cavities 
are indeed required in our models. Thus, our protrusions or their clumps along the SE-NW 
bipolar lobe likely trace the unreso lved dense wal l materials (i.e., swept-up material) of that 
bipolar lobe as suggested in |SC06| and ISSMC07I . They only appears at low velocity likely 
because the SE-NW bipolar lobe is lying in the plane of the sky, as supported by the similar 
brightness of the SE and NW components of the bipolar lobe, in both the HST optical 
and the NIR images. They are faint here in our CO observation probably because their 
CO emission merges with that of the extended circumstellar envelope and halo, and is thus 
mostly resolved out. Two other protrusions are also seen at low velocity along the E-W 
bipolar lobe (Fig. Wfa), tracing the unresolved dense wall materials of the bipolar lobe. 
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5.2. The shaping mechanism and its consequences 



The PPN 19475+3119 is not spherically symmetric, with a torus-like envelope at the 
common waist of the two bipolar lobes. The torus-like envelope is expanding at a speed of ~ 
12.5 km s _1 , and thus has a dynamical age of ~ 1900 yr, with a radius of ~ 5000 AU (~ 1"). 
The E-W bipolar lobe has a mean inclination of ~ 20°, so the mean deprojected velocity is ~ 
27/ sin(20) ~ 80 km s _1 at the tips and the mean deprojected distance of the tips is ~ 21000 
AU (~ 4"3). Thus, the dynamical age of the E-W bipolar lobe is estimated to be ~ 1200 
years (±400 years), more than half of that of the torus-like envelope. The SE-NW bipolar 
lobe is about half the length of the E-W bipolar lobe and may have a similar dynamical age 
to the E-W bipolar lobe within a factor of two. The bipolar lobes and torus-like envelope are 
believed to be physically related. The envelope is believed to be the AGB envelope formed by 
the late AGB wind. The two bipolar lobes could trace two bipolar outflows produ ced by two 
distinct b ipolar post- AGB winds ejected in two different directions, as proposed in lSCOa and 
SSMC071 . The densed-wall structures with tenuous interiors (cavities) of the bipolar lobes 
could result from the interact ions between the post- AGB winds and the e nvelope as shown in 
hydrodynamical simulations ( jLee fc Sahaill2003l ; iLee. Hsu, fc Sahaill2009l ). Also, according to 
our models, the envelope, originally either a torus or a spherical shell, is disrupted by the two 
bipolar lobes or outflows, and it appears as it is now because of the disruptive interactions. 

Could the torus-like envelo pe originally be a torus, like the de nse equatorial waist struc- 
ture often seen in other PPNe flVolk et al.ll2007l ; ISahai et al.ll2008l )? If so, it could be due to 
an eq uatorially enhanced mass loss in the late AGB phase (see the review by lBalick fc Frank 
20021 ). Such a dense waist, where large grains can grow, seems to b e needed to account for 
the millimeter /submillimeter excess toward this source, as argued by lSSMC071 . This toroidal 
envelope may help collimate the E-W bipolar outflow in its pole direction. However, it is 
hard for the toroidal envelope to collimate two bipolar outflows that are almost perpendicu- 
lar to each other, such as those in this PPN. Instead, the envelope appears to be disrupted 
by the two outflows, as discussed above. 

On the other hand, could the torus-like envelope be a remnant of a spherical AGB 
shell in the inner part ? The tor us-like envelope has an expansion velocity similar to that of 
the surrounding halo (1SSMC07I ). and thus may have the same origin as the halo, which is 
the spherical circumstellar shell formed by an isotropic AGB wind in the past. As discussed 
before, a spherical shell model with 2 pairs of cavities along the two bipolar lobes can produce 
a torus-like envelope as observed. Thus, the torus- like en velope could indeed be the remnant 
of the AGB shell in the inner part. Note that ISC06I have also tried a similar spherical 
shell model but only with a pair of cavities along the E-W bipolar lobe. As mentioned in 
Section I4.4[ however, this model would produce two clumps that are perpendicular to the 
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E-W bipolar lobe, but not at the common waist of the two bipolar lobes as observed. A 
similar model with 2 pairs of cavi ties has been adopted to exp lain the clumpy hollow shell 
morphology of the Egg Nebula by iDinh- V-Trung fc Liml ( 20091 ). In that PPN, the spherical 
AGB envelope is also disrupted by the outflows. In this PPN, the spherical shell is disrupted 
more dramatically, leaving behind a torus-like envelope. 

The envelope, originally either a torus or a spherical shell, will not be able to collimate 
the two bipolar outflows, as discussed. Therefore, the two collimated outflows in this PPN 
must be produced by two post- AGB winds that are intrinsi c ally collimated or jets, as pro 



posed in other PPNe and young PNe flSahai fc Traugerlll998l : ISahaill200ll ; iLee fc Sahaill2003 



One popular model to produce a collimated jet requires a binary comp anion to accrete 
the m a terial from the AGB star and then launch the jet (see the review by iBalick &: Frank 
20021 ). iGledhill et al.l (120011 ) also suggested that a binary interaction is needed to produce 
the two-armed spiral structure seen in the NIR polarization image of this PPN. The binary 
interaction can also transform the isotropic mass loss of the AGB wind into an equatorial 
enhanced mass loss (IMastrodemos &: Morrislll998l ). producing a toroidal envelope as needed 
in our model. However, it is unclear h ow two jets can be launched in this PPN to produce 
the two bipolar outflows (see ISSMC071 . for various possibilities). 



6. Conclusions 



IRAS 19475+3119 is a quadrupolar PPN, with two bipolar lobes, one in the east-west 
(E-W) direction and one in the southeast-northwest (SE-NW) direction. The E-W bipolar 
lobe is known to trace a bipolar outflow and it is detected at high velocity. The SE-NW 
bipolar lobe appears at low velocity, and could trace a bipolar outflow moving in the plane 
of the sky. Two compact clumps are seen at low velocity around the common waist of the 
two bipolar lobes, spatially coincident with the two emission peaks in the NIR, tracing dense 
envelope material. They are found to trace the two limb-brightened edges of a torus-like 
circumstellar envelope, expanding away at ~ 12.5 km s _1 . This torus-like envelope can be 
reproduced reasonably well with either a torus or a spherical shell, both with two pairs of 
cavities along the two bipolar lobes. Here, the torus could come from an equatorial enhanced 
mass loss and the spherical shell could come from an isotropic mass loss, both in the late 
AGB phase. In either case, the circumstellar envelope appears to be disrupted by the two 
bipolar outflows. 



We thank the referee for the valuable comments. M.-C. Hsu appreciates valuable discus- 
sion with Chun-Hui Yang, especially in the development of the radiative transfer code. C.-F. 
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Table 1. The best fit parameters for our models 



Parameters 


Torus 


Torus 


Torus 


Shell Model 




no cavities 


+ EW cavities 


+ 2 pairs of cavities 


+ 2 pairs of cavities 


V 


8 


8 


8 





Inclination % 


15 


15 


15 




equator P.A. 


15 


15 


10 




n (cm -3 ) 


7.5 x 10 3 


7.5 x 10 3 


9.0 x 10 3 


4.3 x 10 3 


r (K) 


28 


28 


28 


28 


v e (km s _1 ) 


12.5 


12.5 


12.5 


12.5 


M (M yr^ 1 ) 


1.9 x 1(T 5 


1.9 x 1(T 5 


2.3 x 10~ 5 


2.7 x 10" 5 



Note. — The torus model with p=8 has a half opening angle of 9q ~ 24°. Here the 2 pairs 
of cavities include the E-W pair and SE-NW pair of cavities. The E-W pair of cavities have 
a P. A. =88° and a mean inclination angle of ~ 20°, with the western cavity titled toward us. 
The SE-NW pair of cavities have a P.A. =—40° and are assumed to have an inclination of 0°. 
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Fig. 1.— CO J=2-l contours of IRAS 19 475+3119 superposed on top of (a) the HST 
image obtained with F606W filter ( 1SSMC071 ) and (b) the NIR image (with yellow contours, 



Gledhill et al.ll200l[ ). The central source position is marked by a "+". The synthesized 
beam of the CO observation is shown in the left corner, with a size of l"16xl"01. The 
CO intensity maps integrated over three selected LSR velocity intervals show the torus-like 
structure (black contours) , redshifted (red contours) and blueshifted (blue contours) outflow 
lobes. Contour levels all start from 2a with a spacing of la, where a is 0.45 Jy beam -1 km 
s -1 . Labels A-E show the cuts for the position- velocity (PV) diagrams in Figure [3j 
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Fig. 2. — Channel maps of CO J=2-l emission in IRAS 19475+3119. The velocity width is 
3.0 km s -1 . The systemic velocity is 18.25±1 km s _1 , as found in our model. The contour 
levels start from 2a with a spacing of la, where a =0.06 Jy beam -1 . The stellar position 
is marked by a cross. The synthesized beam is shown in the first channel map at the lower 
right corner. The velocity in each channel is shown in the top left corner. 
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Fig. 3. — PV diagrams of the CO emission, with the cut directions shown in Figured^. The 
contour levels start from 2a with a spacing of 2a, where a = 0.06 Jy beam -1 in (a)-(d) and 
0.03 Jy beam -1 in (e). The spatial and velocity resolutions are shown as a rectangle in the 
bottom right corner, (a) Cut A (along P.A.=15°) goes through the emission peaks of the 
two low-velocity clumps, showing a ring-like PV structure, (b) Cut B (along P.A. =— 10°) 
goes along the axis half way to the edges of the clumps. It also shows a ring-like structure 
but with a weaker emission around the systemic velocity, (c) Cut C (along P. A. =—30°) goes 
through the edges of the clumps, showing a breakup of the ring-like PV structure, (d) Cut 
D (along P.A.=— 75°) goes perpendicular to the axis connecting the two clumps, showing 
the emission shifted to the east on the blueshifted side and to the west on the redshifted 
side, (e) Cut E (along P. A. =88°) goes along the high- velocity component in the east- west 
direction, showing a velocity increasing with the distance from the source for the outflow, 
and a hint of bifurcation of the PV structure (solid lines). 
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Fig. 4. — CO spectrum toward the center position, averaged over a region with a diameter 
of 2"5. Note that a bigger region will include emission from the two bipolar lobes. The 
vertical dashed line shows the systemic velocity. The spectrum shows more emission on the 
redshifted side than the blueshifted side, as expected for an expanding envelope. 
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Fig. 5. — Comparison of our torus model with the observed torus-like envelope in (a) inte- 
grated intensity map as shown in Figured! (b) spectrum as shown in Figure HJ and (c)-(f) 
PV diagrams as shown in Figure |3k-d. Black contours and spectrum are from the model. 
The best fit model requires p = 8 and an inclination of ~15°. 
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Same as Figure [5] but for the torus + EW pair of cavities model. 



- 22 - 



< 



CL 



CJl 
CD 




Aa (*') 



J 


V 


A. 
J 


b) 


_l T -I j 

"V i 



-20 -15 -10 -5 



10 



15 20 



QJ 
(f) 

6 

c 
o 

'w 

o 

LL 




South 

-i — i — | — i — ^ 

North 




East 



-12 



-4 



Velocity relative to the systemic velocity (km/s) 
Fig. 7. — Same as Figure [5] but for the torus + 2 pairs of cavities model 
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Fig. 8. — Same as Figure [5] but for the spherical shell + 2 pairs of cavities model. 



